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We propose a theoretically consistent modification of gravity in the infrared, 
which is compatible with all current experimental observations. This is an 
analog of the Higgs mechanism in general relativity, and can be thought of 
as arising from ghost condensation — a background where a scalar field has 
a constant velocity, (0) = M^. The ghost condensate is a new kind of fluid 
' that can fill the universe, which has the same equation of state, p = —p, 

as a cosmological constant, and can hence drive de Sitter expansion of the 
^ ■ universe. However, unlike a cosmological constant, it is a physical fluid 

d ' with a physical scalar excitation, which can be described by a systematic 

effective field theory at low energies. The excitation has an unusual low- 
energy dispersion relation u"^ ~ /M"^. If coupled to matter directly, it 
gives rise to small Lorentz- violating effects and a new long-range spin 
dependent force. In the ghost condensate, the energy that gravitates is not 
the same as the particle physics energy, leading to the possibility of both 
sources that can gravitate and antz-gravitate. The Newtonian potential is 
modified with an oscillatory behavior starting at the distance scale Mpi/M^ 
and the time scale Mpj/M^. This theory opens up a number of new avenues 
for attacking cosmological problems, including inflation, dark matter and 
dark energy. 



1 Introduction and Summary 



Gravity remains the most enigmatic of all the fundamental interactions in nature. 
At the Planck scale, large quantum fluctuations signal a breakdown of the effective 
field theory description of gravity in terms of general relativity, leading to the ex- 
pectation that gravity and our notions of space and time must be radically altered 
at short distances. Meanwhile, continuing experimental probes of gravity at large 
distances reveal strange phenomena at many length scales, from the flattening of 
galactic rotation curves to the accelerating universe. Traditionally, these phenomena 
are explained by invoking new sources of matter and energy, such as dark matter 
and a tiny cosmological constant. It is however worth contemplating the possibility 
that gravity itself is changing in the infrared in some way that might address these 
mysteries. This approach is further motivated by the cosmological constant problem, 
which seems to be associated with extreme infrared physics. 

As a first step in this program, it is an interesting theoretical challenge to find 
a theoretically consistent and experimentally viable modification of gravity in the 
infrared (IR). Such attempts have been made in the past PHEllllElinilZlElinilinillll- 

Many of these theories are in one way or another varieties of scalar-tensor theories 
of gravity, and therefore are not true modifications of gravity. For example, the 
the long-range interactions between two masses in these theories is at most changed 
by a constant factor. Much more interesting ideas are those that modify gravity 
interactions in a dramatic way at large distances, such as the Fierz-Pauli theory of 
massive gravity and the Dvali-Gabadadze-Porrati (DGP) brane- world models with 
infinite volume extra dimensions However, it is not yet entirely clear that these 
theories are experimentally viable. In both of these theories, there is an additional 
scalar degree of freedom, vr, interacting with matter with gravitational strength. In 
the case of massive gravity it is the helicity zero longitudinal mode of the graviton 
[121 Cni mj; and in the DGP model it is the "brane-bending" mode [T3] (see also 
[in]). Because of the presence of the vr's, the theory reduces to scalar-tensor theory of 
gravity, rather than general relativity, at distances smaller than the length scale of IR 
modification. Furthermore, the vr's become strongly coupled in the ultraviolet (UV) 
at a scale Auv determined in terms of the IR modification scale Air and the Planck 
scale Mpi. At Auv, the effective theory breaks down and the physics is sensitive 
to the unknown UV completion of the theory. That such a UV scale should exist 
is not surprising. After all, even in non-Abelian gauge theories, if we modify the 
theory in the IR by giving the gauge boson a mass m, this introduces a new degree 
of freedom, the longitudinal component of the gauge field. This degree of freedom 
becomes strongly coupled in the UV at scale of order iirm/g, determined by the 
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scale of IR modification and the gauge coupling g. However, the situation is both 
qualitatively and quantitatively different in the gravitational case. The highest strong 
coupling scale Auv can be pushed in all of the above cases is [TH HH] , 

Auv ~ (A?rMpi)1/1 (1.1) 

For Air ~ Hq, today's Hubble scale, we have Ayy ~ 1000 km, which is much larger 
than the distances to which we have probed gravity. Roughly speaking, this happens 
because gravity is being coupled to a sector that is sick in the limit as gravity is 
turned off (Mpi —>■ oo). The vr's have no kinetic term in this limit, even though they 
have cubic and higher order self-interactions. The theory only becomes tenuously 
healthy due to the coupling with gravity: mixing with gravity generates small kinetic 
terms for n of the correct sign (in flat space), but leads to strong coupling physics at 
low energies. This does not necessarily mean that massive gravity or the DGP model 
cannot describe the real world, only that it appears necessary to make nontrivial 
assumptions about the UV completion of the theory. 

In this paper, we present a fully controllable and calculable theory where gravity 
is modifled in the infrared, in a way that avoids these strong coupling issues. It is 
useful, though not necessary, to view this as arising from "ghost condensation." A 
ghost condensate is a new kind of fluid that can flll the universe. It has the equation 
of state p = —p, just like a cosmological constant A, and can therefore gives rise to 
a de Sitter phase of the universe even if A = 0. But the ghost condensate is not a 
cosmological constant, it is a physical fluid with a physical scalar excitation. This 
background can be thought of as arising from a theory where a real scalar fleld is 
changing with a constant velocity 

(0) = or {(j)) = M\ (1.2) 

and there is a physical scalar excitation tt around this background, 

^ = MH + 7i. (1.3) 

We assume that the fleld (p has a shift symmetry 

0^0 + a, (1.4) 

so it is derivatively coupled, and the quadratic lagrangian for (p begins with a wrong- 
sign kinetic term [in (H ) signature] 

C^ = -^d''<pd^(l) + --- . (1.5) 
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The wrong-sign kinetic term means that the usual background with (0) = is un- 
stable. However, as is familiar from tachyon condensation, we will see that higher 
order terms can stabilize a background where 7^ 0. Note that this background 
breaks Lorentz invariance spontaneously — there is a preferred frame where is spa- 
tially isotropic. This is not fundamentally different from the way in which the cosmic 
microwave background radiation (CMBR), or any other cosmological fluid, breaks 
Lorentz invariance. What is novel about the ghost condensate is that, unlike other 
cosmological fluids (but like a cosmological constant), it does not dilute as the universe 
expands.^ 

We will construct a systematic derivative expansion governing the low-energy 
physics of the excitation n about the background. As we will see, these excitations 
have two important properties, that combine to give rise to a number of surprising 
new phenomena in the IR. 

First, the low-energy effective action for vr has the universal form 



5~ d^x 



^*^-^(VM= + -]. (1.6) 

SO that the tt's have a low-energy dispersion relation of the form 

(Throughout this paper k represents the spatial 3-momentum and V represents the 
spatial gradient. In this introductory section, we will not be careful about 0{1) 
coefficients in front of various terms.) Ordinarily the dispersion relation for light 
excitations is of the form tu^ oc k'^. Here, there is no k"^ term at all, so the excitations 
do not move with any fixed speed. Instead, the group velocity is v"^ ~ k'^/M'^, so 
that larger lumps of vr move more slowly than smaller lumps. Note that this is not 
a Lorentz-invariant dispersion relation, as is to be expected since the background 
breaks Lorentz invariance. 

Second, there are two distinct kinds of energy in these backgrounds: a "parti- 
cle physics" energy and a "gravitational" energy, and they are not the same. This 
happens because the background (0) = M^t breaks time translational invariance. It 
also breaks the shift symmetry, but the diagonal combination is left unbroken to 
generate the "time" translations in this background. The Noether charge associated 
with this unbroken symmetry is what an experimentalist coupled to tt would call the 



^Attempts to modify gravity with a dynamical preferred frame violating Lorentz invariance have 
been made, e.g., most recently in Ref. j^, but with different physical conclusions. 
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Fig. 1. Excitations n about the ghost condensate have a dispersion relation cu'^ oc k'^, 
so that small lumps of tt expand more quickly than larger lumps of tt. 



conserved energy in this background, and this "particle physics" energy density takes 
the form 

and is manifestly positive for small fluctuations. On the other hand, the energy 
density that couples to gravity is Tqo, which is associated with the broken time- 
translation symmetry, and is not the same as £pp. Indeed, 

^grav = 7})0 ~ M^TT + • • • (1.9) 

begins at linear order in n, therefore lumps of n can either gravitate or an^i-gravitatc 
depending on the sign of tt! In fact, the "particle physics" energy is nothing other 
than the inertial mass, while the "gravitational energy" is the gravitational mass, 
so we can say that the tt excitations maximally violate the equivalence principle. It 
is worth mentioning that excitations of tt with exotic gravitational properties are 
not easily produced if tt couples only to gravity. In fact, we find that turning on tt 
gravitationally requires long distances and time scales, as discussed below. 

The standard model fields do not need to couple directly to the ghost sector. How- 
ever, if there is a direct coupling, then the leading interaction (a derivative coupling 
of (p to the axial vector currents) gives rise to potentially observable consequences, 
including a splitting between particle and antiparticle dispersion relations, and a new 
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Fig. 2. Lumps of tt can either gravitationally attract test particles (in regions where 
TT > 0) or repel them (where tt < 0). 



spin-dependent inverse- square law forces mediated by vr exchange. This inverse square 
law behavior is a direct consequence of the uj'^ oc k"^ dispersion relation. 

Even if the standard model fields have no direct couplings to the ghost sector, 
they will indirectly interact with it through gravity, and the propagation of gravity 
through the ghost condensate gives rise to a fascinating modification of gravity in the 
IR. As an example, consider a cloud of dust that quickly collapses to form a star or 
a planet at some time t — 0. At linear order, the gravitational potential around this 
planet is modified at a length scale Tc, but only after waiting a much longer time tc, 
where 

Mpi Ml 

Starting at time tc the potential becomes oscillatory at distances of order Tc, and the 
modification reaches out to much larger distances r ':$> rc only after waiting an even 
longer time t {r/rc)tc. Furthermore, in a flat background, the ghost condensate has 
a Jeans-like instability at long wavelengths, so that the modulation of the potential 
begins growing exponentially at late time t ^ tc- As with ordinary fluids, however, 
this Jeans instability disappears due to Hubble friction in a cosmological background 
(such as de Sitter space) with H > t~^. In this case, at late time t 3> tc, the potential 
looks hke an oscillatory modulation of 1/r such as shown in Fig. 3. 

There are a number of interesting ranges of parameters for the mass scale M 
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Fig. 3. Example of a late-time Newtonian potential 

associated with the ghost condensation. If the ghost condensate is to give the observed 
acceleration of the universe today, even with A = 0, then naturalness suggests 

M ~ IQ-^eV. (1.11) 

In this case, the distance scale where gravity is modified is ~ -f^o^^, but the time it 
takes to see this modification is tc ^ Hq^\ So no modifications of gravity can be seen 
directly, and no cosmological experiment can distinguish the ghost-driven acceleration 
from a cosmological constant. If the standard model has direct coupling to the ghost 
sector, we can still hope to detect tiny Lorentz-violating effects and long-range spin- 
dependent inverse square force laws. It is also possible that M is much larger; this 
could be the remnant of ghost condensation triggering an earlier de Sitter epoch in 
the universe, to drive inflation jTTj. For instance, we can have the parameters 

M~10MeV, tc^H^^ 1000 km, (1.12) 

so that we begin to start seeing an oscillatory modulation of Newtonian gravity at a 
distance of ~ 1000 km, but only recently in the history of the universe! 

This physics is an exact analog of the Higgs mechanism for gravity. Recall that 
general relativity is a gauging of spacetime translations. In particular the time dif- 
feomorphism symmetry is spontaneously broken by our background (0) = M'^t, and 
therefore we should expect gravity to become "massive." There are two differences 
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with the famihar gauge theory case. First, our background breaks Lorentz invariance, 
and therefore we do not get the massive graviton of a Lorentz invariant theory, which 
would have 5 polarizations. In fact we are only adding one additional physical degree 
of freedom, which can be thought of as the scalar excitations it. Second, while in 
the gauge theory case we get a Yukawa-like exponential suppression of the potential 
at large distances, in our case we get an oscillatory modulation instead — basically 
because this is a negative mass-squared term for the graviton. This difference in sign 
is to be expected in all comparisons of gravity and gauge theory, physically arising 
for the same reason that like charges repel in gauge theory and attract in gravity. 
It is also analogous to the fact that fluctuations of a charged fluid acquire a posi- 
tive plasma frequency Up > 0, while gravity induces an < signalling the Jeans 
instability. 

We have described the physics so far in the language of gravity propagating 
through a new fluid. This is like the description of massive gauge theories in terms 
of propagation through the Higgs condensate. As familiar from the gauge theory, 
there is also a "unitary gauge" description of the physics, where the excitation vr is 
set to zero ("eaten"), and where the interpretation of the theory directly as an IR 
modification of gravity is clear. The unitary gauge is arrived at by choosing the time 
coordinate t such that 

^{t,x) = MH. (1.13) 

That is, is taken to be the clock. This is equivalent to making a diffeomorphism to 
set TT = 0. In this gauge, there are terms in the action that are not fully diffeomor- 
phism invariant. At quadratic level for the fluctuations h^^, = g^y — r]^^ around flat 
space these new terms in the lagrangian take the form 

AC = aoM'^hlo - aiM'^KijKij - a2M'^{Kiif, (1.14) 

where ^ 

Kij = 2^dohij - dihoj - djhoi), (1.15) 

and ao,i,2 are 0{1) coefficients. The fully non-linear action will be presented in section 
6. These terms are clearly "mass" terms and unusual kinetic terms for h^^,. 

It is also instructive to understand the counting of degrees of freedom in this the- 
ory in the (3 + l)-split ADM language ^H]; where the metric is split into N ^f^", 
Ni ^ goi, and 'jij ^ gij. Usually in GR, both N and Ni are not dynamical fields, they 
are Lagrange multipliers, and their equations of motion give 1 + 3 = 4 constraints on 
the , reducing the 6 degrees of freedom in •jij to the two physical polarizations of the 
graviton. Our background breaks time-diffeomorphisms, but leaves (time-dependent) 
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spatial diffeomorphisms intact. As a consequence, the Ni are still Lagrange multi- 
pliers enforcing 3 constraints. However, N is no longer a Lagrange multiplier, and 
the N equation of motion then simply determines N. There are therefore 6 — 3 = 3 
physical degrees of freedom — we have added one degree of freedom. Like all gauge 
"symmetries," diffeomorphisms are not a symmetry, they are a useful redundancy of 
description. The physical issue is what the degrees of freedom and their dynamics 
are. The additional degree of freedom is not manifest in the unitary gauge descrip- 
tion of the physics, and it is most convenient to explicitly introduce the Goldstone 
field TT that non-linearly realizes the symmetry. In this way, arriving at this theory 
through the background (0) = JVPt of a scalar is the analog of a linear sigma model 
description of gauge symmetry breaking. This is not necessary however — we could 
have simply started with the theory of the tt and its (non-linear) transformations 
under time diffeomorphisms explicitly, just as we can directly write the chiral lagran- 
gian for pions without any reference to an underlying linear sigma model. This line 
of thought suggests another way in which our model might arise from a more funda- 
mental underlying theory. One of the well-known conceptual difficulties in putting 
quantum mechanics and gravity is the problem of time: quantum mechanics picks 
out a direction of time, in conflict with general covariance. Obviously this problem 
must be resolved at some zeroth order, however, if any remnant of this clash remains 
at low energies, it will be described by our theory at long distances and times. 

In the rest of the paper we will discuss this physics in detail. In section two 
we study ghost condensation in general, first in the absence of gravity, and then 
in a cosmological Priedmann-Robertson- Walker (FRW) background. As we will see 
Hubble friction in an expanding universe is crucial in determining the background. 
We discuss how de Sitter phases of the universe can arise from ghost condensation in 
section 3. In section 4, we discuss the low-energy effective action for the fluctuations 
around the ghost condensate in the absence of gravity, and perform the correct power- 
counting analysis to show that the theory is (non-trivially) a sensible, quantum- 
mechanically stable effective theory at low energies, so the tt's are a healthy sector 
even in the absence of gravity. In section 5, we discuss the signals that would arise 
from a direct coupling of the standard model fields to the ghost sector. In section 6, 
we systematically construct the low-energy effective theory for the ghost condensate 
coupled to gravity. In section 7, we study the modification of gravity in the ghost 
condensate in flat space at linear order, arising from mixing between the scalar sector 
of gravity and the tt excitations. Here we uncover the Jeans instability of the fluid 
at large wavelengths, and the oscillatory modulation of the Newtonian potential at 
large distances and late times. We also make some brief comments on the interesting 
question of corrections to non-linear gravity. In section 8, we repeat this analysis 
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for a de Sitter background, where the Jeans instabihty is removed. We end with a 
discussion of some of the many open questions and avenues for further research. 



2 Ghost Condensation 



2.1 Generalities 



Consider a scalar field with a shift symmetry (j) t—>- (j) -\~ a ensuring that it always 
appears derivatively in the action. Assume that it has a wrong-sign quadratic kinetic 
term, i.e., it is a ghost field: 



Actually, the overall sign of the lagrangian is completely unphysical, so what is ac- 
tually relevant is that once a choice has been made for the sign of the lagrangian 
for, say, the standard model and gravity sector, a ghost is a field with opposite sign 
kinetic term. The difficulty is then that, with both kinds of fields, the energy of the 
theory is unbounded from below, and therefore the vacuum {(p) = is unstable. Vac- 
uum decay will inevitably happen quantum-mechanically. Consider the minimal case 
where (p is only coupled gravitationally. Kinematically, nothing prevents the vacuum 
from decaying to a pair of 0's and a pair of gravitons. There is more and more phase 
space for this for higher energy gravitons (and more negative energy 0's), but if the 
theory has a UV cutoff at a scale M, then the rate for vacuum decay per unit volume 
is r ~ Afjy/MpY- For Auv ^ 10"'^ eV, this decay will happen within our Hubble 
volume within a Hubble time. Other possible decays and constraints were considered 



Let us compare this with the more familiar example of a tachyon, which is a scalar 
with negative mass squared. In this case the vacuum (0) = is also unstable. Even if 
classically we can sit at the point = forever, quantum-mechanically we inevitably 
roll off the top of the potential with a time scale of order m~^. In the tachyon case, we 
are accustomed to the idea of tachyon condensation: if there are higher order terms 
that turn the potential around, e.g., V{(j)) = — |m^0^ + |A0^ + ■ ■ ■ with A > 0, then 
there can be a stable vacuum where (0) 7^ 0. 

We can ask whether it is possible to have find a stable "ghost condensate" vacuum 
for a theory with wrong-sign kinetic terms. In fact, general relativity can be viewed as 
such a theory. The Vierbein and spin connection form the gauge fields for IS0{3, 1), 
the group of Lorentz transformations and translations on 4-dimensional Minkowski 
space. The gauge group is noncompact, so fiuctuations about the vacuum = 
include ghosts, but fiuctuations about nontrivial vacua such as fiat space (e*^) = 5" 



ldyd^<j) + ■■■. 




in Ref. ^U^. 
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Fig. 4. The function P. The shaded region is the one for which the small fluctuations 
about the background have healthy two-derivative time and spatial kinetic terms. 

are stable. In our theory with a single scalar, we can imagine that there are other 
terms such as (90)^ in the lagrangian that "stabilize" the ghost in the same way 
the (f)^ term stabilizes the tachyon. One immediate difference with the tachyon case 
is that all of these higher order terms are non-renormalizable; if these terms are to 
become as important as the (90)^ term as they must in order to stabilize the theory, 
then it appears that a description in terms of effective field theory is not possible. As 
we will see this is in fact not the case, and there is a systematic low-energy expansion 
for fluctuations around the ghost condensate. But let us press on for a bit, and start 
by assuming that the lagrangian has the form 

C^M^P{X), (2.2) 

where X = d^(j)d^^(j). We have now chosen (j) to have the dimension of length, so that 
X is dimcnsionless, while M is a constant with the dimension of mass. Note that we 
are neglecting terms involving more than one derivative acting on 0; we will come 
back to these terms below. A typical form for the function P we are considering is 
shown in Fig. 4. 

The equations of motion that follow from this lagrangian are 

d,[P'{X)d^(t>]-^- (2.3) 
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From this we see that we obtain a solution provided that = constant. If d^cf) is 
time-hke, there is a Lorentz frame where 

= ct, (2.4) 

where c is a dimensionless constant. These solutions are actually the maximally 
symmetric solutions, since they break the (f) shift symmetry and the time translation 
symmetry down to an unbroken diagonal shift symmetry. 

Let us now consider small fluctuations about these solutions 

(p = ct + Tx. (2.5) 

The lagrangian for quadratic fluctuations is 

C = M^[ [P'{^) + 2^P"{^)\ TT^ - P\c^) (V7r)2} . (2.6) 

We see that the fluctuations have the time and spatial kinetic terms of the usual signs 
provided that c is such that 

P'{c^) > 0, P'{c^) + 2^P"{^) > 0. (2.7) 

Around these backgrounds the excitations of vr have positive energy, and so any weak 
coupling of this sector to other sectors (such as gravity and the standard model) will 
not suffer from the violent UV instabihty associated with the ghosts around the usual 
background with 9^0 = 0. This does not uniquely fix c, any value of c satisfying these 
conditions is stable against small fluctuations (see Fig. 4). 

More generally, we have additional terms in the effective lagrangian arising from 
terms with two or more derivatives acting on 0: 

C = M^P{X) + M^Si{X){U(t)f + M^S2{X)d^d''(l)d^d,(t) + ■■■ . (2.8) 

(We have assumed for simplicity here that there is a t— > —0 symmetry that forbids 
terms with an odd number of 0's in the lagrangian). With these additional terms 
included, there is still a solution for any value of c, since the equation of motion 
always takes the form 



function of dch, d"^ 



0. (2.9) 



The quadratic lagrangian for vr now includes terms with more derivatives on the vr's. 
We will see below that Hubble friction drives us to a background where the spatial 
kinetic terms coming from P' vanish and that the leading spatial kinetic terms are 
fourth order and arise from the S\^2 terms in Eq. ()2.8j) . 
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We can now understand why it is possible to have a good low-energy effective field 
theory for the physics, despite the appearance of the non-renormalizable interactions 
P{d4>y, S{{d(f)Y) and so on. The point is that the background sits at some specific 

value of c, and therefore does not sample large regions of the functions P, 5'i,2, 

Terms with n > 1 derivatives acting on give terms of order 9"7r, which are higher 
order in the derivative expansion. After Taylor expanding the coefficient functions 
5'i,2j ■ ■ ■, we obtain an effective lagrangian with a finite number of parameters that 
controls small fiuctuations about a given background. 

Another way to look at this is to note that the background is characterized by 
a symmetry breaking pattern. For example, the ghost condensate breaks Lorcntz 
transformations down to spatial rotations. The tt field is the Goldstone boson and 
must non- linearly realize the full Poincare symmetry of the theory, and the tt effective 
lagrangian is highly constrained by symmetry considerations. This point of view will 
be more fully developed below. 



2.2 Hubble friction 

So far, we have a continuum of possible backgrounds labelled by the parameter c. 
The situation changes in the presence of gravity and an expanding universe. Gravity 
is coupled to the effective lagrangian of the previous section via the usual minimal 
coupling: 

C = y/^M^P{X) + ■■■, (2.10) 

where X — g^"dn(f)dy(f). Again, we begin by ignoring the higher order terms 5'i^2, ete. 
The energy-momentum tensor of the ghost condensate is 

T^, = -M^P{X)g^, + 2M^P\X)d^ct>d,(t>. (2.11) 

In general, we expect the background to be an expanding FRW universe, with metric 

ds"^ = de - a\t)dfl^, (2.12) 

where dfP is the spatial metric for a maximally symmetric 3-dimensional space. For 
(f) depending only on t, the (j) equation of motion in these coordinates is 



d. 



a^P'{4>'^)(f> =0. (2.13) 



Note that because this is the equation of motion, it holds even if some other sector 
is dominating the expansion of the universe. 
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From here we see that 



0P(02) ^ ^0 as a ^ oo. (2.14) 

a'^{t) 

So, in the far future as a ^ oo, we are driven to either 

= or P'(02) = 0, (2.15) 

depending on the initial condition. The first option is what we are famiUar with from 
massless scalar fields in an expanding universe — the scalar velocity rcdshifts to zero 
as the universe expands. However, the point with = is the one corresponding to 
the unstable theory with the wrong sign kinetic term. The other possibility is that 
the theory is driven to <j) = c*, with P'{cl) = 0. If P"{cl) > 0, the fluctuations about 
this point are stable with suitable higher derivative terms as we will see momentarily. 

These conclusions are robust under inclusion of the higher-order terms 5*1^2, 

In the case where the universe becomes flat in the asymptotic future, (p — c implies 
that V ijVv<t> — 0, and so all of the higher order terms vanish in the equation of 
motion and we are still driven to the point where P' = 0. If the asymptotic metric 
is de Sitter space with a Hubble constant iJ, it is no longer true that for <f> = c, all 
higher covariant derivatives vanish; for instance, = 3Hc, so the new terms make 
non- vanishing contributions to the equation of motion. In this case, we are not driven 
to P' = exactly, but the important point is that Hubble friction singles out some 
value of c* with P'(c^) ~ powers of H/M. These are small for H/M <^ 1, which is 
required by validity of the low-energy effective theory. 



2.3 Low energy dispersion relation 

We have seen that gravity forces us into an extremum P'{cl) = (for an asymptoti- 
cally flat spacetime). Recall that the spatial kinetic terms for tt from the expansion 
of P about the background arc proportional to P', so if P' = there is no (V7r)2 tcr m 
in the quadratic lagrangian for tt. However, the fluctuations can still be stabilized by 
(V^tt)^ terms, which arise from the higher order terms proportional to 5'i,2, 



Si{cl)+S2{cl) (VV + ---- (2-16) 



This stabilizes fluctuations about (f) — c^t provided that Si{cl) + S2{cl) < 0. The 
quadratic lagrangian for tt is then 

jC = - ^M'(V'7^)^ (2.17) 
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where we have put P"{cl) = j, which can always be achieved by rescahng between 
P and M^, and we have defined = -2M'^{Si{c^)^ + ^2(0^))- The tt modes then 
have the low-energy dispersion relation 

= jr.k\ (2.18) 

This conclusion also holds for an asymptotically de Sitter space. As an example, 
consider the lagrangian of the form 

C = ,/^M^[P{X) + Q{X)R{n<l))]. (2.19) 

For depending only on t, the equation of motion ()2.13p is modified to 

dt {a^ [(P'(0') + Q'(0')i?(0 + 3/70))20 - dt{Q{^^)R'{4> + 3if0))] } = 0. (2.20) 

As a 00 we have 

^'(0^) + Q'{<P'^)R{4> + SH<i))] 20 ^ (2.21) 

Q{(j)'^)R'{4> + 3H(p) const. (2.22) 

Expanding in powers of vr, we find that the coefficient of the (Vvr)^ term is exactly 
—M^{P'{(j)'^)+Q'[(j)'^)R[(j)+3H(f))), which is driven to zero (for 7^ 0) by the expansion 
of the universe. The generality of the fc^ dispersion relation is easier to see from the 
non-linear realization formulation, to be discussed in section 6. 

If the universe did not start in the special configuration <p = c*, the (Vtt)^ term 
will not vanish exactly today and its size depends on the history of the universe. 
As we will see later, even in the absence of the (Vtt)^ term there will be a k"^ piece 
induced by mixing with gravity in the dispersion relation. The coefficient of the (Vtt)^ 
term can be made negligibly small compared with the gravity-induced piece if it was 
driven to zero e.g. by a period of infiation in the early universe. On the other hand, 
it would also be interesting to study the effects of a significant (Vvr)^ term that may 
still remain in the current universe. 

2.4 Two kinds of energy and (anti)gravity 

Our background 

(0) = c,t (2.23) 

spontaneously breaks time translation invariance. It also breaks the shift symmetry 
— i> + c. However, a diagonal sum of these two shift symmetries is left unbroken 
as the time translation in the background. The combination 

£pp= (d^xT^o-c^Q (2.24) 
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of the Noether charge Jd^xToo associated with time translations, and the charge Q 
associated with the shift symmetry, is the conserved energy in this background. 
The usefulness of Spp is that it is the conserved energy that an experimentalist would 
infer from the dynamics of the tt fluctuations 



(2.25) 



It is the conservation of this quantity that ensures the stability of the vacuum against 
small fluctuations. On the other hand, this is not the energy Tqo that couples to 
gravity. Indeed, from the expression for T^^ we can infer that 

Too = c^M^ + • • • . (2.26) 

This begins at linear order in tt. Therefore, Tqo can be either positive or negative, 
depending on the sign of tt, and locahzed excitations of the ghost condensate can give 
rise to antigravity! Together with the unusual dispersion relation for tt, this will be 
the root of the exotic infrared phenomena we will find below. 



3 de Sitter expansion from ghost condensation 

Our background spontaneously breaks Lorcntz invariance, since there is a preferred 
frame in which is spatially isotropic. This is not conceptually any different from the 
Lorentz violation due to any cosmological fluid that fills the universe, like the cosmic 
microwave background radiation. Indeed, because of Hubble friction it is natural 
that this frame coincides with the frame in which the cosmic microwave background 
is isotropic. There is however an important difference between the ghost condensate 
and conventional forms of matter and radiation: it does not redshift away as the 
universe expands. Indeed, in the far future the energy-momentum tensor becomes 

T^. ^ -g^,.M^P{cl), (3.1) 

i.e., precisely the form of a cosmological constant! If -P(c^) < 0, the universe will 
approach a dc Sitter metric. However the ghost condensate is not a cosmological 
constant, it is a physical fluid with a physical scalar excitation. 

The cosmological constant in this theory is simply A = — M^P(0)/M|j. Even if 
this vanishes, the ghost condensate can drive the current acceleration of the universe 
if M^P(c^) ~ — (10~^ eV)"^, and assuming that there are no very large or small 
dimensionless parameters in c, this implies that M ~ 10"'^ eV. It is also possible 
for the ghost condensate to drive an early inflationary phase of the universe, as will 
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be discussed in Ref. jTJj. The fact that accelerated expansion of the universe can be 
driven by the kinetic terms was discussed before, under the names of "k-inflation" |2I] 
and "k-essence" |^. In those works, is required to move over a large range of the 
kinetic function, which is truncated to low orders or assumed to be of some special 
form. These models therefore require some knowledge of the fundamental theory and 
are not based on consistent effective field theories. 

If is not exactly at the solution = where P'{c\) = 0, then we have 

P'{^') ~ 4- (3-2) 

From the energy- momentum tensor of the ghost condensate Eq. (j2.1H) . we see that 
a nonzero P' couples to gravity like non-relativistic matter, namely Tqo ~ M'^P', 
Tij = 0. Therefore in the approach to its final position at = c*, the ghost condensate 
sources gravity like non-relativistic matter. This raises the prospect that the ghost 
condensate may contribute to dark matter as well as the vacuum energy, hinting that 
these phenomena may both have their origin in a modification of gravity at large 
distance scales. This possibility will be investigated in future work. 



4 Low-energy effective theory 



We would now like to examine the low-energy effective theory for vr in more detail, in 
the absence of gravity. We will consider the theory expanded around the background 
P'{cl) = 0, which we are driven to by Hubble friction. We show that this theory is 
truly healthy in the absence of gravity, in contrast to the analogous scalar sectors in 
DGP and massive gravity, which only become healthy due to the coupling to gravity. 
In particular, we will address the regime of validity of the effective theory and the 
question of stability. 

Expanding X = cl + 2c^i: + — (Vvr)^, we can see that the low-energy effective 
action for the tt's begins as 

S= Id'xdt ^M^TT^ -^M'^iV^ny (4.1) 



(4.2) 



where we have suppressed the higher order terms that arise from the 5*1^2 terms giving 
rise to the (V^vr)^ terms. 

Clearly for this to be a sensible low-energy effective theory, the cutoff has to 
be beneath the scales M, M, which we will take to be roughly comparable. If we 
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canonically normalize the tt^ kinetic term, we find that the (V^tt)^ term and the 
cubic interaction 7r(V7r)^ term are both mass dimension 6 operators. Which term is 
more important at low energies? A closely related question is that of stability. The 
(particle physics) energy is independent of the 7r(V7r)^ term because it is linear in tt, 
but the quartic term (Vtt)^ contributes a negative-definite term to the energy density 
that will be larger than the positive contribution from the (V^vr)^ term provided that 
Vtt > k/M. It is therefore possible to find directions in field configuration space with 
arbitrarily long wavelengths in which the energy can be made more and more negative. 
This means that we have to address the nonlinear stability of the theory. This is most 
pressing at quantum level: what stops quantum fluctuations from lowering the energy 
and leading to vacuum instability? 

The key to understanding all of these issues is to do the proper power-counting 
analysis for the scaling dimension of all the operators in our effective theory (see 
Ref. for a very clear discussion of these issues). In a relativistic theory, the 
scaling and mass dimension of couplings are the same, while this is not the case in 
the non-relativistic theory we are studying. Let us begin then by identifying the 
scaling dimensions as the energy is scaled hy E sE, where s is some constant, or 
in terms of time t ^ s~^t. We determine the scaling dimension of space x and the 
field TT by requiring that the quadratic part of the action is invariant. This then fixes 

E ^ sE, (4.3) 

t ^ s-H, (4.4) 

X (-^ s^^^'^x, (4.5) 

TT s^/V. (4.6) 

The scaling of x is expected since u (x k"^ , as to ^-^ scu we should take k — > s^^'^k and 
so X 1-^ s^^/^x. With this scaling, it is easy to check that the leading interaction 

d^xdiM^niVny (4.7) 

scales as s^^^. It is therefore (just barely!) an irrelevant operator, becoming less and 
less important at low energies compared to M. All the other operators are even more 
irrelevant than this leading one. This implies that there is a regime of low energies 
and momenta, and low field amplitudes, where the expansion is under control.^ In 
particular, the quantum fluctuations are not large enough to make the higher order 

^This is a non-trivial result. For instance, suppose it had turned out that the leading spatial 
kinetic term for tt was of the form (VV^tt)^, then the scaling dimensions of the field would have 
been different {x h^- s~^^^x and tt s°7r), and the leading operator 7r(V7r)^ would have scaled as 
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terms compete with the quadratic terms. There are therefore no large quantum 
instabihties in the IR; only in the UV where the theory becomes strongly coupled is 
there a danger of instabilities from the higher order terms. This is however sensitive 
to the UV completion, and in this respect the situation is similar to QED or any other 
theory that gets strongly coupled in the UV. In any case, if we work with a UV cutoff 
A that is somewhat smaller than the scales M, M, then the strongest dimensionless 
coupling strength in the theory is set by powers of A = A/M. The timescale for any 
quantum instability, if it exists, will be exponentially small in 1/A, and can easily be 
longer than the age of the universe. 

For completeness note also that there is an operator that is not forbidden by any 
symmetries, namely (Vtt)^. This is definitely a relevant operator, however, as we have 
seen, Hubble friction always drives the background to a point where the coefficient of 
this operator vanishes. 

The fact that the theory is healthy independently of gravity is in contrast with 
case of the massive gravity or the DGP model. In our model there is a limit where 
the scalar decouples, while keeping the strength of gravity fixed. To do this, we take 

Mpi oo and also all matter sources T^,y oo, keeping the gravitational scale 
Tfj^iy/Mpi fixed, while also keeping M fixed. In massive gravity or DGP model, there 
is a scalar coupled to ordinary matter with the same strength as the graviton and its 
effects do not vanish in this limit. 

The scales M that we will be considering will turn out to be quite a bit smaller 
than the TeV scale, in the range 10^^ eV < M < 10 MeV, and the n sector will 
need to be embedded in a UV completion above this scale. However, gravity and the 
standard model are very weakly coupled to the tt sector, so the breakdown of the tt 
effective theory at M will not appreciably affect the standard model sector. 

5 Direct coupling to Standard Model fields 

The standard model fields need not have any direct couplings to the ghost sector 
(other than indirectly through gravity). However, it is interesting to consider possible 
direct couplings of the standard model fields to 0, preserving the shift symmetry on 
(p. The cutoff of the standard model sector can be much larger than the cutoff of the 
ghost sector as long as the couplings between them are weak enough. The leading 

s~^/^ and been relevant, becoming strong and leading to a breakdown on the eflFective theory at low 
energies. Fortunately this did not happen for us. It is amusing to note that our leading operator is 
more irrelevant for all spacetime dimensions > 4, while for d = 3 it is classically marginal! In this 
case the 1-loop beta function must be computed to determine the IR behavior of the theory. 
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derivative coupling of to standard model fermions is familiar from Goldstone boson 
interactions, and is of the form 

^C = Y.''-^ija^^d,(t^. (5.1) 

(Here we have reverted to canonical normalization for the cf) field.) A systematic 
discussion of all possible Lorentz violating effects for standard model fields in flat 
space can be found in Ref. [211 123 • Note that these couplings can be removed by field 
redefinitions on the fermions '0 — *■ e^'^^f'^^^^-ip; however, if the symmetry ip e«c^^'^ 
is broken by mass terms or other couplings in the lagrangian, the interaction cannot 
be removed. If ipjip'^ are paired with a non-zero Dirac mass term moipip'^, then the 
vector couplings (with + c^c = 0) can be removed, while the axial couplings 

AC ~ 4^^'7V*<9^0 (5.2) 

r 

remain. Expanding (p = M'^t + tt, we have 

1 - 

AC ~ /i^'7V* + t^^t'^T^^^m^' where fJ^ = (5.3) 

The first term is a Lorentz and CPT violating term, that gives rise to a different 
dispersion relation for particles and their antiparticles. In particular, a left helicity 
particle and right helicity antiparticle have dispersion relations of the form 



uj = ^J{\p\±iiY + ml (5.4) 

where the + sign is for the left-helicity particle and the — sign is for the right- 
helicity antiparticle, and the signs are reversed for a right-helicity particle and left- 
helicity antiparticle [IHIEII. Furthermore, if the earth is moving with respect to the 
background in which the condensate is spatially isotropic, there is also an induced 
Lorentz and CPT-violating mass term of the form 

/^^77^^ . vearth- (5.5) 
In the non-relativistic limit, this gives rise to an interaction Hamiltonian 

IJ'S ■ {/earth- (5.6) 

The experimental limit on /i for the coupling to the electrons is of the order < 10~^^ 
GeV |2Hj, and to the proton and the neutron are < 10"^"^ GeV [211 EHI, assuming 
l'5'earthl ~ 10"'^. A review of Lorentz and CPT experimental tests and extensive 
references can be found in Ref. jHlj . 
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Also strikingly, the exchange of tt gives rise to a new long-range force. In the 
non-relativistic limit, we have a derivative coupling to spin 

ACr^l;S-VTT, (5.7) 
r 

As with usual Goldstone bosons, the exchange of vr produces a long-range spin- 
dependent potential. Ordinarily for Goldstone bosons this produces a poten- 
tial, but here, because of the dispersion relation instead of the usual k"^ dispersion 
relation, we get a 1/r potential: 

^ M^F^ r ■ ^ ■ ^ 

Therefore these goldstones mediate spin-dependent inverse-square law forces! Note 
that we have imagined exactly static sources in deriving this force law, ignoring the 
retardation effects — the u'^ piece in the denominator of the vr propagator. For massless 
particle exchange in relativistic theories, retardation effects can be ignored as long 
as the sources are static on a timescale longer than the time it takes light to travel 
between them. In our case, however, since u"^ oc k^, for sources a distance r apart the 
static limit force law we have found is only valid on timescales longer than 



'^r^Mr\ (5.9) 



We have focused here on direct couplings to the (p field that are linear in (j). 
These can all clearly be forbidden by a — —0 symmetry. There are however some 
couplings between the standard model fields and that are inevitably generated by 
graviton loops. The leading operators of this type have the form 

^ 0^'^d,4>d,^, (5.10) 



where is some dimension 4 standard model operator. Putting in the background 
= then gives rise to tiny dimensionless Lorentz violating effects of size ~ 
M^/Mpi, for instance differing speeds of light for different particle species. Also, for 
^SM ~ T^^'^, this operator can act as a coherent source for 0, however because of the 
Mpi suppression, this is never competitive with the existing terms in the lagrangian 
for 0. 

There is a large literature on Lorentz violation in the standard model [221 • The 
rotationally invariant effects have been thought of as arising from a non-zero value 
for the time component of a background vector field A^, (Aq) ^ 0. Our model is 
of this type with ^4^ — > 9^0, and therefore provides a fully consistent framework to 
study Lorentz violation, including gravitational effects. 
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6 Ghost condensate coupled to gravity 

In this section, we construct the effective theory of the ghost condensate coupled to 
gravity We can simply do this by minimally coupling our lagrangian for to gravity, 
but it also instructive to see how this works directly in a "unitary gauge," where the 
modification of gravity is seen transparently 

To define unitary gauge, note that we can use general coordinate invariance to 
use as a time coordinate, i.e., we choose coordinates {t,x^) such that 

(f){t,x)^t. (6.1) 

This is a good gauge choice for small fluctuations about a time-dependent solution 
such as the one we are considering. This choice of gauge eliminates as a degree of 
freedom, and leaves residual gauge freedom corresponding to time-dependent spatial 
diffeomorphisms : 

t^t, x'^x''{t,x). (6.2) 

We can write the unitary gauge lagrangian in terms of quantities that transform 
simply under these residual diffeomorphism symmetry. 

Let us first see this works at quadratic level in the action expanding around flat 
space, Qfj^i, — r]fjti,-\-hni,. Under general diffeomorphisms generated by x^ — > x^-\-^^{x), 
we have as usual 

5K, = -{d^i, + d,Q. (6.3) 
Under the residual unbroken diffeomorphisms generated by ^i, we then have 

5hQQ = 0, 5hQi = -9o6, 5hij = -{di^j + dj^i). (6.4) 

What invariants can we use to build the action? Recall that we are assuming that 
flat space is a good background here, so any action must begin at quadratic order in 
the /i's. Clearly, since /ioo is invariant, the leading term is of the form 

Jd^xdt^M'^hlo- (6.5) 

As we have mentioned, because these terms do not preserve the full diffeomorphism 
invariance of the theory, we are really introducing an additional degree of freedom. 
To see this explicitly, it is convenient to re-introduce the field tt that restores the full 
diffeomorphism symmetry of the theory — we can achieve this simply by performing 
a broken ^° diffeomorphism and promoting ^° = tt to a field. Then 

hoo hoo - 2(?o7r, hoi hoi - diir, hij hij, (6.6) 
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and we see that the M*^/iqq action has generated a time kinetic term iM*^7r^ but no 
spatial kinetic term for vr. In fact, it is easy to see that only adding the M^/iqq term 
to the action, we have not modified GR at all (at least classically), since this is a 
partial gauge-fixing term that fixes to /iqo = gauge. We must therefore go to higher 
order, to make vr a dynamical field, and see a real modification of GR rather than 
just a gauge fixing.^ 

In this language, it is easy to see why we do not get a (Vvr)^ kinetic term. We could 
get this from terms of the form h^^ in the action, since upon introducing the tt we have 
hoi — > hoi — diH. However, K^^ is not invariant under the residual diffeomorphisms! 
There are invariants that can be constructed, but they involve higher derivatives, for 
instance 

= ^i^ohij - djhoi - dihoj) 

^{dohij - djhoi - dihoj + didjir) (6.7) 

is invariant (and is the linearized extrinsic curvature of constant time surfaces in the 
theory). The leading terms at quadratic level in the effective theory are of the form 



S= d^xdt 



C2 2 



(6.8) 



Terms linear in K can be forbidden by time-reversal invariance. If we do add for 
example hooK -k'V'^iT + ■ ■ ■, it contributes a term uk"^ to the vr dispersion relation. 
However, this does not qualitatively change the physics as ~ fc^ still holds. Terms 
linear in tt do not contribute to the Hamiltonian and hence do not affect stability 
either. 

Upon re-introducing tt, ()6.8|) gives us the spatial kinetic terms, 

1 



5- = d-^xdt 



M' = {M' + M"), (6.9) 



so we see that the k^ spatial kinetic terms are a direct consequence of the assumption 
that flat space is a good background, and the unbroken residual spatial diffeomor- 
phisms of the theory. 

We can repeat the same analysis for a de Sitter background, the metric is of the 
form 

ds^ = (1 + hoo)df - e^^^dx'dx^{6ij - hij) + 2hoidtdx\ (6.10) 

^The Hqq term does not fix time reparameterization globally. In particular, static solutions are not 
unique. The correct procedure is to define static solutions as the limit of time-dependent solutions, 
as will be done in specific examples below. 
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and the linearized diffeomorphism transformations arc 

Shoo = -25oCo, Shoi = -di^o - e^^^do^i, Sh^ = 2H^oSij - di^j - dj^i, (6.11) 

where = Vtii^C"- As before, we can restore full diffeomorphism invariance by pro- 
moting ^0 to a dynamical field tt. Then the only invariants under linearized diffeo- 
morphisms are 

hoo - 27r, and - e'^^^dihoj + djho, - didjir) + 26ijH'k. (6.12) 

As before, we cannot write an invariant term that contains a (Vtt)^ kinetic term. 

Let us now proceed to the systematic construction of the effective theory at non- 
linear level. First, consider the scalar quantity 

X = g>'''d^(l>d,(t>. (6.13) 

In unitary gauge, 

X ^ (6.14) 

and we see that in unitary gauge the lagrangian can be an arbitrary function of g^^. 
Writing the lagrangian in terms of X is useful because it is obvious how to generalize 
to an arbitrary gauge. 

We can also define the unit vector perpendicular to the surfaces of constant 0: 

This can be used to project any tensor into tensors with all indices parallel to the 
surfaces of constant (j). For example, the induced metric is simply 

l^iv = -{Qhv - n^n^), (6.16) 

which satisfies n'*7^,^ = 0, and is therefore a metric on the surfaces of constant 0. 
We can use this to write other tensors such as the intrinsic curvature R{'j) and the 
extrinsic curvature 

Kf^u = 7/V,n„ (6.17) 

which is a symmetric tensor satisfying n'^K^i^ — 0. 

Since the residual diffcomorphisms depend on time, do is not a covariant derivative. 
The covariant time derivative on tensors with all indices parallel to surfaces of constant 
4> is 

L>'^T-... = ^(n^V^r-...)ll, (6.18) 

V A 
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where '||' denotes the projection onto surfaces of constant using ra^. (The normal- 
ization is dictated by the fact that n^d^ = in unitary gauge.) 

The most general effective lagrangian invariant under spatial diffeomorphisms in 
unitary gauge is therefore 

£ = v^[f(X) + lQi{X)K^ + \Q2{X)K'^K,, + • • • 

(6.19) 

+ terms involving time derivatives . 

Here the ij indices are raised and lowered with , 7*-' . Again we have assumed the 
time-reversal symmetry forbidding the terms linear in K. Note also that at this order 
we could have written down a term involving the intrinsic curvature R^'^\ however, 
the full four-dimensional curvature i?*^^-* is given by i?*^^) = R^^^ + Kf^ — K^, so we 
can always write a term linear in R^^^ as a usual Einstein-Hilbert term plus terms 
proportional to K'^, Kfj. 

Note that we can easily restore the dependence on the ghost fluctuations vr by 
using the general formulas above. Using the relation ^/—g = 1/7/ VX, we see that 

F{X) = M^^^. (6.20) 



Let us now consider the expansion about fiat space in unitary gauge. The terms 
that do not involve derivatives in the unitary gauge lagrangian are 

C = . 6.21 



We would like to have a solution where g^y = r]^^, so X = r^^^ = 1. It is easy to see 
that this requires -P(l) = to avoid tadpoles for fluctuations of the spatial metric, 
and -P'(l) = to avoid tadpoles for g^'^. The fact that P{1) = can be interpreted 
as the condition that there is no cosmological constant (see Eq. (12.111) ). The fact 
that P'{1) = implies that the ghost fluctuations have no /c^ term in their dispersion 
relation as we have seen, and the leading terms are the ones (see Eq. (j2.6|l ). 

Summarizing, the leading low-energy effective action is, in unitary gauge, 

"1 1 ~ 1 ~ 



S = Jd^xdt^ ^M^{X - If - -M^K^ - -M'^KijK'^ + 



(6.22) 



7 Infrared modification of gravity 



We now consider the modification of gravity in the ghost condensate, working at 
linearized level. At this order, the only effect comes from the mixing between tt 
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and gravity, and so it suffices to only consider scalar perturbations of the metric. A 
general scalar perturbation can be parameterized by 



/ioo = 2< hoi^diX, hij ^ 25ij<lf' + didjTi. (7.1) 

Using the the diffeomorphisms generated by — ^diT), Co — X ~ |^o?7, we can bring 
the metric to the longitudinal gauge, 

/loo = 2$, hoi = 0, hij = 2^dij. (7.2) 

We can analyze the physics at the level of the quadratic action (or equivalently 
the linearized equation of motion). Let us begin by working at the level of the action. 
The Einstein-Hilbert action in this gauge is, 

-3(^0^)^ - 2V* • V$ + (V^)^! , M|i = (7.3) 



Note that despite the appearance of the time derivative on ^ in this action, neither $ 
nor ^ are dynamical fields. Indeed it is easy to go to momentum space and diagonalize 
the 2x2 kinetic matrix mixing $, \E'; the zeros of the determinant of this matrix give 
the dispersion relation which is just k'^ = 0. This is in keeping with the fact that the 
only propagating fields in gravity are the two polarizations of the spin-two massless 
gravitons. Note also that in the Newtonian hmit, o;^ <^ A;^, so the action reduces to 
-2V* • V$ + (V*)^, and the equation of motion then fixes = 

We now include the ghost condensate. We again work in the Newtonian limit 
where co'^ <^ k'^. Then at quadratic level the lagrangian for tt becomes 

jC^^M^^-tt)^ --M\V^7r)\ (7.4) 

We can then put ^ to its equation of motion = and go to canonical normahzation 
$ = $c/(v^Mpi), TT = TTc/M^, to find the lagrangian 

We can see that there are small mixing terms between $c and tTc- Going to momentum 
space, the lagrangian in the (tTc, $c) basis is 



where we have defined 



M2 2 

m = —;= , a = — r. 

V2Mpi M2 



(7.7) 
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We see that the vr dispersion relation is modified from lo^ = a^k'^/M'^ due to 
mixing with gravity. The dispersion relation is found by setting the determinant of 
the 2x2 matrix above to zero, which gives 



2 2 2 2 2 , 2 



uj' = - = - ^^k\ (7.8) 

M2 M2 M2 2M|i ^ ' 

For very low momenta k < m, we find c<j^ < 0, signaling an instability. This is the 
analog for our fiuid of the usual Jeans instability for ordinary fiuids; for a fiuid of 
mean density p and pressure p, the usual Jeans instability shows up as a modification 
of the dispersion of the form^ 

= ~ ^^^^^ = ^i^- ^^■^'> 
In our case, the largest imaginary magnitude of uj for the instability is 



M3 



^inst = = (7.10) 

Just as the usual Jeans instability is removed by Hubble friction in an expanding 
universe, we will see in the next section that for our case that in a de Sitter background 
with H > r, this instability is also removed. 

As we will now see, F"^ is also the timescale over which modifications of gravity 
take place, at a length scale of order m~^. To see this, we can look at the modification 
of the $c propagator. If we immediately go to the static limit — > 0, we see that 
the $c propagator is simply modified from 1/fc^ — * l/(^^ ~ '"^^) (remembering that 
k is the 3-momentum). This looks like a negative mass squared for $c (the sign is 
dictated by the requirement of a healthy tt time kinetic term), and so gives rise to an 
oscillatory sin(mr)/r or cos(mr)/r modification of the Newtonian potential, rather 
than the familiar e~"^^/r Yukawa form. However, as we already saw for the spin- 
dependent force mediated by vr exchange, because the tt's do not move at the speed 
of light but instead slow down in their propagation at larger length scales, retardation 
effects are very important, and we cannot immediately go to the static uj limit. It 
is straightforward to find the ($c'^'c) propagator by inverting the 2x2 kinetic matrix, 
yielding 



^ a^ni?k'^ 



M'^uj^ — a'^k^ + a'^m?k'^ 



x(-^)- (7-11) 



^Strictly speaking, the dispersion relation Eq. 1)7. 9|l is not correct, because in ordinary gravity 
flat space is not a solution when p 7^ (the 'Jeans swindle'). In the ghost condensate, flat space is 
a solution, so the dispersion relation Eq. (|7.8|) is not a swindle. 
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The factor in the brackets represents a change in the Newtonian potential. To get an 
(9(1) modification of gravity, we must have 

This tells us that we have to be close to "on-shell" for the tt excitations; directly in 
terms of lu we then have 



UJ — 



M 



So we see that the width of the region in u; is small, set by Alo ~ F, and therefore we 
must wait for a long timescale tc ~ F""^ to see any modifications of gravity. It is also 
easy to see that the natural length scale of the modification is Tc ~ m~^. Indeed, if 
we work with dimcnsionlcss variables ili, k defined by a; = 2Fa), k — mk, we see that 
the modulation factor for the Newtonian potential is simply 



1- 



Cj'^-k^ + k^ 



(7.14) 



and so modifications occur for u), A; ~ 0(1). Wc could have also seen that r^tc are 
the natural length and time scales for modification also directly from the action. By 
further defining $ = m$c and tt = F^^tTc, the action can be seen to be 0(1) in all 
hatted variables. 

We therefore conclude that modifications of gravity take place at a distance Tc 
and time tc given by 

Mpi . Ml, 

The fact that tc ^ makes sense from a number of points of view. Note that 
the scale Vc is completely determined by the M^/iqq part of the lagrangian. However, 
as we have already remarked, with only this term, classically we have not modified 
gravity at all, only fixed to hoo = gauge. Only with the other terms (which lead 
to the k^ spatial kinetic terms for tt) do we see that there is a genuine modification 
of gravity — to see any change in the potential we then have to wait until both terms 
become important, and the timescale is then larger than what one might expect from 
the coefficient of the h^Q term alone. 

Let us see how the modification works explicitly in an example, by calculating 
the effective gravitational potential that is felt by a test mass, outside a source which 
turns on at the origin at time t — 0, i.e., a source 

Pm{r,t) = S^{r)e{t). (7.16) 
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Of course energy is really conserved here; this can be thought of as an approximation 
to a situation where a cloud of dust quickly collapses to a point-like mass at the 
origin. We are interested in the gravitational force that a test particle feels a distance 
r away, as a function of time t > 0. Given the $c propagator we have calculated, it 
is easy to Fourier transform back to position space to find 

*(r,t) = -^[l + /(r,t)], (7.17) 

where I{r, t) is a spatial Fourier integral over the momentum k. Introducing the 
dimensionless variables 



we have 



/(r,t) = -J dup^^(l-cos\i{TuVl^) 

-K \ Jq [U-^ — U) ^ 

+ r du^^^ (l - cosiTuV^^f^)) \ . (7.19) 

Jl [U'^ — u) ^ ^1 

For times t <ti tc{T <til), I{r, t) — > and there is no modification of the standard 
Newtonian potential. However, for i > tc{T > 1), there are modifications. The 

presence of the exponentially growing cosh(T-uVl — u'^) term in the first integrand is 
a refiection of the Jeans-type instability we found in fiat space. For large T, therefore, 
the first integrand will dominate the integral for I{r,t). The exponential in the cosh 
is maximized for u = 1/a/2; approximating the integrand by a Gaussian around this 
point and performing the resulting Gaussian integration then gives us an excellent 
approximation for I{r,t): 



/(M)-^exp^-- + -jsm^^j. (7.20) 

So, ioi R <^ T {r I r c ^ 1 1 tc) I we indeed find an oscillatory behavior for the Newtonian 
potential, which is growing exponentially in time as exp(T'/2) as a consequence of the 
instability! However, for R ^ T, the correction vanishes; it takes a time of order 
t = tc X (r/rc) for any appreciable changes to happen at distance r. The exact 
behavior of this Newtonian potential at various time is shown in Fig. 5-7. 

Note that for M ~ 10^'^ eV, which could drive the acceleration of the universe 
today, Tc ~ Hq^ is about the size of the universe, while tc ^ -f^o"^ is much larger than 
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Fig. 6. Same as the previous figure except for t — (2r) ^. Note the axes have been 
rescaled. 
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Fig. 7. Same as the previous figure this time for t = 5r~^. Note again that the 
axes are rescaled and the correction at hnear level has become (exponentially) larger 
in amplitude. This is due to the Jeans instability of the fluid in flat space 

the age of the universe! Thus, modifications of the Newtonian potential will occur 
only at much later time than the present age of the universe. On the other hand, for 
M ~ 10 MeV, we can have tc ~ Hq^ comparable to the present age of the universe, 
and these oscillatory modifications of gravity take place on a distance scale of ~ 1000 
km! 

We close this section by making some brief comments on the modifications of 
gravity at non-linear level. In backgrounds with large gravitational fields such as a 
neutron star, large redshift factors will make vary by 0(1), and we might wonder 
whether this leads to large deviations from GR. Recall that there is a decoupling limit 
in which all the effects of the ghost condensate on the matter sector vanish: Mpi and 
all matter source energy densities pm ^ oo, while keeping the usual gravitational 
scales pm/Mp^ fixed, and also M fixed. In this limit the tt's have no back-reaction 
at all on the metric. Beyond this, it is easy to see that with only the P{X) term in 
the lagrangian, if we begin with initial data with P'{X — cl) — (as we know we 
are driven to by Hubble friction), then X — cl for all time. Thus, in a non-trivial 
gravitational potential, (j) adjusts itself so that X = g'^''dfj_(j)du(f) is fixed and there is 
no modification of GR. Only with the additional terms in the action (that give rise to 
the A;^ spatial kinetic terms) do we get a deviation from X = c^. As we saw already 
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at linearized level, however, this effect takes time to build up. At any rate, due to 
the existence of the healthy decoupling limit, there is some limit on M which will be 
in agreement with all tests of GR. Clearly M ~ 10"^ eV will be fine, since in this 
case even the most naive estimate of the scale of modification is M^/Mpi ~ Hq. A 
more detailed investigation is needed to find the precise bounds near the much larger 
values of M ~ 10 MeV, that lead to deviations from the linearized Newtonian gravity 
at timescales of order the current age of the universe at distances of order ~ 1000 km. 
It is possible that nonlinear solutions evolve to states outside the regime of validity 
of the effective theory, but the same arguments lead us to expect that this takes very 
long times. 



In the previous section we studied the ghost condensation and modification of gravity 
in a Minkowski background, and we saw that there is a Jeans-like instability. As in 
the case of the Jeans instability for ordinary matter, we expect that it will go away in 
an expanding universe if the expansion rate exceeds the growth rate of the instability. 
In this section, we study the ghost condensation and modification of gravity in such 
an expanding de Sitter background. This is relevant since the current expansion of 
the universe is accelerating and we are approaching a de Sitter space. One could also 
consider the possibihty that ghost condensation is responsible for inflation in the early 
universe. In that case, the fluctuations of ghost condensation in de Sitter space can 
be relevant for generating the density perturbations that seed structure formation in 
the universe. 

It is more convenient to do the analysis using the linearized equations of motion, 
so for ease of comparison with our results from the last section, let us begin by writing 
the linearized equations of motion in flat space. In longitudinal gauge, the Einstein 
tensor in flat space is 



8 The Ghost Condensate in de Sitter Space 



00 



(8.1) 






and the equations of motion are 



Goo + 2 ^'(^00 - 29o7r) = 0, 
-Goi - a^M\diK) = 0, 



(8.4) 
(8.5) 
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Gi, + a^M^doK6ii = 0, 



1 a2y2 

-9o(9o7r - -hoo^ 



K = 0, 

2 . 



(8.6) 
(8.7) 
1 for conve- 



where we have kept only h^Q and i^T^ terms and set M|j = 1/ (87rG') 
nience. 

From the third equation, for i ^ j, we find that $ = Further dropping 
subleading terms for M/Mpi <^ 1, we find 

2V2$ + M*^ - M^doTT = 0, (8.8) 
-2(9o<9i$ - a^M'^diV^TT = 0, 
29^$ + a^M^doV^TT = 0, 



- ^0$ + 



M2 



-TT = 0. 



(8.9) 
(8.10) 

f8.11) 



The second and third equations in the above both imply that 

1 



^0$ 



Inserting this into the vr equation of motion yields 



(8.12) 



(8.13) 



From here we can see the modified dispersion relation for vr, and together with the 
equation for V^$, the modified propagator for $ can be derived. 

We repeat this exercise in a de Sitter background. As before, we can use diffeo- 
morphisms to bring the metric to the following form, 

ds"^ = (1 + 2<l>)rft2 - a(t)2(l - 2^)dx^ a{t) = e^* . 



The equations of motion ()8.8p - ()8.1H) are modified to 

-GHdo"^ + 2—$ - 6/7^$ + M\<l> - don) = 0, 
-di {2H^ + 29o^) - a^M^di Udo^ + ^vr ) = 0, 



(8.14) 

(8.15) 
(8.16) 



2Hdo{^ + 3^) + 2(9^^ + 6H\^ + ^) + ^($ - ^) 



■^($-*)-<5, 



a 



' y2 N 

39o^ + — TT 



0. 



(8.17) 
f8.18) 
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We consider that the Hubble expansion parameter hes in the interesting range, 



r = — ^ <H<^m[ = < M. 8.19 

For H <^T the space-time is effectively ffat. On the other hand, if H > m, the length 
scale of the potential modulation is outside the horizon and the timescale for vr to 
modulate gravity is way longer than H~^, so the interesting modification of gravity 
could not have happened. 

Again, from the third equation with i 7^ j, we have $ = Eliminating $ and \I/ 
from the above equations and neglecting higher order terms, we obtain the equation 
of motion for the vr field. 

By Fourier transforming the spatial dependence, —k"^, the equation of each 

mode k can be written as 

d',n + H i^^^-f^] don-^^ U - ^ V = 0. (8.21) 

The Sovr term provides the friction (anti-friction) for modes with k/a>m or k/a< 
m /V3 (m/v^ < k/a < m). The redshift will eventually bring each mode to the 
regime k/a <^ m/^/S and the friction term becomes simply HdoTi. For H > T, 
the friction term dominates and the instability disappears just like the usual Jeans 
instability. 

To study the modification of gravity due to the ghost condensate, we should 
examine the gravitational potential by eliminating vr instead. The equation for the 
gravitational potential $ without additional source is 

(a^ + 3H8, + 2H^) * + ^ % + 5* = 0. (8,22) 

or for each momentum mode k, 

(9o^ + 3Hd, + 2H^) $ + ^ f ^) $ - ^ f ^) $ = 0. (8.23) 

They are equivalent to Eqs. ()8.20p . ()8.21|) . Absence of instability for i7 > F can also 
be shown with the $ equation. At a large distance (^ m~^) from a local lump of 
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excitations, one can neglect the last two terms of the equation. The solutions at a 
fixed comoving distance r are simply 

$(r, t) ~ bi{r) e"^* + b2{r) e"'^*. (8.24) 

They decay at least as fast as the redshifting so there is no growing potential. 

The equation of the gravitational potential with ordinary matter source can be 
similarly derived. It is convenient to decompose $ into two parts as 

$ = $GR + $mod, (8.25) 

where each part satisfies 



72 

o^"^"" 2 

and 



^$GR = ^, Peff = p-4V2p, (8.26) 



4r^ /V72\2 .p2 y2 



mod 

4r2 v2 

— — $GR, (8.27) 



respectively, for the general scalar-type matter source 

Too = P, To. = dig, a-^T,, = p6i, + (^2ViVj - l^ij"^^) P- (8.28) 
To derive the above equations we have used the conservation equations 

p + 3Hp + 3Hp -g = 0, 

g + 3Hg-^V^p-p = 0. (8.29) 

Note that $gr is exactly the gravitational potential in general relativity and, thus, 
$jnod represents the modulation of gravity due to the ghost condensation. Equation 
fl8.27|) shows that the timescale of the modulation is indeed F"^ for the length scale 

To obtain the potential as a function of the physical distance, we can switch to 
the physical coordinates X = a{t)x, then the equation for $mod becomes 

{do + HX- Vx)'<f mod + 3H{do + Vx)$mod + <l>^od 

4r^ 2 4r'^ 

+ ^(^x) -^-mod + ^ V^$n.od = rPeff (8.30) 
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The modulation of the gravitational potential from a source also starts at the time 
scale as the vr field slowly reacts to the gravitational potential. If we are interested 
in the potential at very late time, we can look for time independent solutions for $ by 
setting 9o = as there is no instability in this case. Assuming spherical symmetry, 
the time independent equation is 



where a = m\X\. 

We first consider the solutions without matter source and, thus, without $gr- The 
regularity at the origin a = requires that all odd-order derivatives of $mod vanish at 
(7 = 0. Hence, there are two independent solutions. The two independent solutions are 
shown in Fig. IHlfor H/T = 2, where we plot the potential multiplied by the distance 
from the origin. One can see that apart from the oscillatory modulation, the potential 
is proportional to at large R. This can be interpreted as arising from an effective 
mass produced by a lump of vr excitations around the origin. The expansion of the 
universe keeps these vr excitations stationary (tt 7^ 0), resulting in a static potential. 
Any linear combination of these two configurations represents a possible gravitational 
potential produced by vr without matter source, and the effective mass can be either 
positive or negative which corresponds to a gravitating or an antigravitating object 
respectively. 

A time-independent solution of the potential with a regular matter source can 
also be obtained similarly. A solution with pes/m^ = 0.1e~^'°'^^ for H/T = 2 is shown 
in Fig. ini Again we see the interesting modulation of the 1/R potential. The mass 
— i?$ [33J seen by an observer far away (m~^ <^ i? -C H~^) is not necessarily the 
same as the mass of the matter source due to the screening or anti-screening of the vr 
field. A general time-independent solution with the same matter source is given by 
the sum of the special solution shown in Fig. El and an arbitrary linear combination 
of the two source-free solutions shown in Fig. |H1 The coefficients of the source-free 
solutions cannot be determined by the time independent equation ()8.31|) but should 
be determined by dynamical considerations. 

This situation is in stark contrast to the corresponding situation in general rela- 
tivity. In general relativity, a spherical symmetric spacetime outside matter source 
is locally static and, thus, locally a Schwarzschild geometry. This fact, known as 
the Birkhoff theorem [33] , makes it possible to characterize the geometry outside the 
source just by a mass parameter. In general relativity, this is true even if the geom- 
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Fig. 8. Numerical solutions without matter source (p = 0) for the boundary con- 
dition ($^od(0) = 0, $:;„d(0) = -0.01), and ($^od(0) = -0.01, ^'l^^^{0) = 0) with 
H/r = 2. 
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Fig. 9. A numerical solution with matter source pcs/m^ — O.le"^''^'^^ for the bound- 
ary condition $mod(0) = *mod(0) = and H/T = 2. 
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etry inside the matter source is completely dynamical. Notice also that even in the 
presence of a usual scalar field, the Birkhoff theorem holds up to linear perturbation 
around Minkowski or de Sitter background, because the scalar field does not couple 
to gravity in the linearized level. (It only appears quadratically or higher order in 
T^u-) On the other hand, in ghost condensation, vr reacts to the gravitational po- 
tential produced by matter source and slowly modulates the geometry in such a way 
that the final stationary configuration depends on the history of the whole system. 
This happens even in the linear perturbation level. What Fig. El shows is just one 
of such possible final configurations. The reason why this remarkable modulation 
happens even in the linear perturbation around a de Sitter background is that we 
have a non-vanishing (p in the background. 

9 Discussion and Outlook 

We have presented a consistent modification of gravity in the infrared. It is useful, 
though not necessary, to view this arising from gravity propagating in a ghost con- 
densate. An immediate application of this model is an alternative to a cosmological 
constant for driving a de Sitter phase in the universe. This can be used to drive the 
acceleration of the universe today, and can also be used to for an early infiationary de 
Sitter phase of the universe, as will be discussed in Ref . [T7] . The physics here is very 
different from standard slow-roll infiation, and the model makes sharp predictions 
about the spectrum of density perturbations that can be excluded or confirmed in 
future experiments jTTj. 

If the standard model fields have direct couplings to the ghost sector, there are 
interesting Lorentz- violating signals and inverse-square law long range forces mediated 
by excitations of the ghost condensate. But most interestingly, gravity is modified in 
the infrared in a remarkable way by the ghost condensate, giving rise to the possibility 
of anti-gravity and an oscillatory modulation of the Newtonian potential at late time 
and large distances. 

Our exploration of these models is still in its infancy. The most pressing question 
is to determine the current experimental limits on the parameters of this theory, given 
the unusual new forces and modifications of gravity it entails. It would also be inter- 
esting to explore further the possibility raised in section 3 that the ghost condensate 
may contribute to both the dark matter and the dark energy of the universe. The 
physics of non-linear gravitational effects, and particularly black holes in these back- 
grounds, is clearly interesting to explore. There are also a number of novel avenues 
to explore for early universe cosmology and models with extra dimensions. It would 
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also be interesting to derive our low-energy effective theory from a more fundamental 
UV complete theory, either arising from ghost condensation or in some other way. 
Finally, can the fact that we have a model where the "energy that gravitates" is not 
the "particle physics energy" help with the cosmological constant problem? 
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